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Abstract 

The thermal relic density of neutralino dark matter in the brane world cosmology is studied. 
Since the expansion law at a high energy regime in the brane world cosmology is modified from the 

Ph! 

one in the standard cosmology, the resultant relic density can be altered. It has been found that, 
if the five dimensional Planck mass M5 is lower than 10 4 TeV, the brane world cosmological effect 
^ ■ is significant at the decoupling time, and the resultant relic density is enhanced. We calculate the 

neutralino relic density in the constrained minimal supersymmetric standard model and show that 
the allowed region is dramatically modified from the one in the standard cosmology and eventually 
disappears as M5 is decreasing. We also find a new lower bound on M5 > 600 TeV based on 
the neutralino dark matter hypothesis, namely, the lower bound in order for the allowed region of 
neutralino dark matter to exist. 
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I. INTRODUCTION 



Recent cosmological observations, especially the Wilkinson Microwave Anisotropy Probe 



(WMAP) satellite have established the ACDM cosmological model with great accuracy, 
and the relic abundance of the cold dark matter is estimated as (in 2 a range) 

VcDMh 2 = 0.1126l°;°m- (1) 

To clarify the identity of a particle as cold dark matter is still an open prime problem 
in cosmology and particle physics. The lightest supersymmetric particle (LSP) is suitable 
for cold dark matter, because they are stable owing to the conservation of R-parity. In 
the minimal supersymmetric standard model (MSSM), the lightest neutralino is typically 
the LSP and the promising candidate for cold dark mater. In light of the WMAP data, 
the parameter space in the constrained MSSM (CMSSM) which allows the neutralino relic 
density suitable for cold dark matter has been recently re-analyzed Q]. It has been shown 
that the resultant allowed region is dramatically reduced due to the great accuracy of the 
WMAP data. 

Note that the thermal relic density of the matter depends on the underlying cosmological 
model as well as its annihilation cross section. If a nonstandard cosmological model is taken 
into account, the resultant relic density of the dark matter can be altered from the one 
in the standard cosmology. A brane world cosmological model which has been intensively 
investigated [3] is a well-known example as such a nonstandard cosmological model. The 
model is a cosmological version of the so-called "RS II" model first proposed by Randall 



and Sundrum 



a, 



where our four-dimensional universe is realized on the "3-brane" located 
at the ultra-violet boundary of a five dimensional Anti-de Sitter spacetime. In this setup, 
the Friedmann equation for a spatially flat spacetime is found to be 



r2 Z*GJ, , p 



where 

p = 96ttGM 5 6 , (3) 

H is the Hubble parameter, p is the energy density of matters, G is Newton's gravitational 
constant with M 5 being the five-dimensional Planck mass, and the four- dimensional cosmo- 
logical constant has been tuned to be almost zero. Here we have omitted the term so-called 



"dark radiation", since it is severely constrained by Big Bang Nucleosynthesis (BBN) j^J. 
The coefficient po is also constrained by the BBN, which is roughly given by p^ 4 > 1 MeV 
(or, equivalently, M 5 > 8.8 TeV) Q]. This is a model-independent cosmological constraint. 
On the other hand, as discussed in the original paper by Randall and Sundrum Q , the pre- 
cision measurements of the gravitational law in sub-millimeter range lead to more stringent 
constraint p^ 4 > 1.3 TeV (or M 5 > 1.1 x 10 8 GeV) through the vanishing cosmological 
constant condition. However note that this constraint, in general, is quite model dependent. 
In fact, if we consider an extension of the model so as to introduce a bulk scalar field, the 
constraint can be moderated because of the change of the vanishing cosmological constant 
condition 6]. Hence, we care about only the BBN constraint on p in this paper. 

Note that the p 2 term in Eq. (J2J) is a new ingredient in the brane world cosmology. Since 
at a high energy regime this term dominates and the universe obeys a nonstandard expansion 
law, some results previously obtained in the standard cosmology can be altered. In fact, 
some interesting consequences in the brane world cosmology have been recently pointed out 
[3,0]. Especially, the thermal relic abundance of dark matter can be considerably enhanced 
compared to that in the standard cosmology Q|. 

In this paper, we investigate the brane world cosmological effect for the relic density of 
the neutralino dark matter in detail by numerical analysis. We will show that the allowed 
region for the neutralino dark matter in the CMSSM is dramatically modified in the brane 
world cosmology. In the next section, we give a brief review of Ref. jzj]. In sec. Ill, we 
present our numerical results for the neutralino relic density in the brane world cosmology. 
We will find some interesting consequences. Sec. IV is devoted to conclusions. 



II. ENHANCEMENT OF RELIC DENSITY IN BRANE WORLD COSMOLOGY 

In this section, we give a brief review on the relic density of dark matter in the brane 
world cosmology with a low five-dimensional Planck mass []]]. 

In the context of the brane world cosmology, we estimate the thermal relic density of a 
dark matter particle by solving the Boltzmann equation 

d/Tt 

_ + 3Hn = - {av )( n 2- n l Q ) i (4) 
with the modified Friedmann equation Eq. (J2J), where n is the actual number density of the 
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dark matter particle, ueq is the equilibrium number density, (av) is the thermal averaged 
product of the annihilation cross section a and the relative velocity v. It is useful to rewrite 
Eq. (0J) into the form, 

dY s 



dx xH 



(™)(Y*-Y£ Q) 
A r £L= {*v)(Y*-Yi Q ), [ ■) ) 



in terms of the number density to entropy ratio Y = n/s and x = m/T, where m is a dark 
matter particle mass, A = 0.26 (g*s/g^ 2 ) Mpm, Mp ~ 1.2 x 10 19 GeV is the Planck mass, 
and x t is defined as 

(6) 



4 _ P 

x t = — 



T=m 



Po 

At the era x -C x t the p 2 term dominates in Eq. (J2J), while the p 2 term becomes negligible 
after x 3> x t and the expansion law in the standard cosmology is realized. Hereafter, we 
call the temperature defined as T t = mx^ 1 (or x t itself) "transition temperature" at which 
the expansion law of the early universe changes from the non-standard one to the standard 
one. Since we are interested in the effect of the p 2 term for the dark matter relic density, 
we consider the case that the decoupling temperature of the dark matter (T^) is higher than 
the transition temperature, namely, x t > x d = m/T d . Using Eqs. © and ©, this condition 
leads to 

/ 2 -1 \ 1/6 

M 5 < ( H± m *^— X --*) (7) 
5 - V 30 96ttG d ) K ' 

^4.6xlQ3 T eVf ™ y/ygoV 73 (g) 
VIOOGeV/ \x d J V ; 

Here we have normalized the decoupling temperature by its typical value in our case esti- 
mated as Xd = m/Td — 20, which is the same scale as in the standard cosmology as is 
shown in the appendix. For M 5 < 10 3 TeV, we can expect significant brane world cosmo- 
logical effects. This, indeed, will be confirmed by numerical calculations in the next section 
(see Fig. 3). 

It is easy to numerically solve the Boltzmann equation Eq. (j3J) with a given (av) and 
Xt- Using appropriate approximations, we can drive analytic formulas for the relic number 
density of dark matter |7J. When we parameterize (av) as (av) = a n x~ n with fixed n = 
0, 1, • • ■, for simplicity, we can obtain simple formulas for the resultant relic densities such 



that 



Y(x — > oo) ~ 0.54—^- for n — 0, 
Acr 



™ 2 



for n = 1, (9) 



A(Ji In x 4 

in the limit Xd <C where Xd is the decoupling temperature. Note that the results are 
characterized by the transition temperature rather than the decoupling temperature. It is 
interesting to compare these results to that in the standard cosmology. Using the well-known 
approximate formulas in the standard cosmology 



Y(x — > oo) ~ - — for n = 0, 
Acr 

2r 2 

for n = l, (10) 

Actx 

we obtain the ratio of the relic energy density of dark matter in the brane world cosmology 
(0(6) ) to the one in the standard cosmology (fi( s )) such that 

fl W ~ 0.54 (-?*-) for n = 0, 

for n — 1, (11) 



ft 



21nxt V^OO. 

where x,i( s ) denotes the decoupling temperature in the standard cosmology. Similar results 
can be obtained for n > 2. Note that the relic energy density in the brane world cosmology 
can be enhanced from the one in the standard cosmology if the transition temperature is 
low enough. With the help of Eqs. (j3J) and (JHJ), we find that the enhanced energy density is 
proportional to M^ 3 ^ 2 (roughly M^ 3 ) for n = (n = 1). 

This is the main point discussed in Ref. Q]. If the above discussion is applied to de- 
tailed analysis of the relic abundance of neutralino dark matter, we can expect a dramatic 
modification of the allowed region in the CMSSM. 



III. NUMERICAL RESULTS FOR THE NEUTRALINO RELIC DENSITY 



In this section, we present the results of our numerical analysis. We calculate the relic den- 
sity of the neutralino Q x h 2 in the CMSSM with the modified Friedmann equation Eq. (J2J). For 
this purpose, we have modified the code DarkSUSY ^ so that the modified Friedmann 
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equation is implemented. In evaluating the relic density, coannihilations with neutralinos, 
charginos, and sfermions are taken into account. 

The mass spectra in the CMSSM are determined by the following input parameters : 

mo, mi/a, A, tan/5, sgn(/x), (12) 

where m is the universal scalar mass, mw2 is the universal gaugino mass, A is the universal 
coefficient of scalar trilinear couplings, tan/5 is the ratio of the vacuum expectation values 
of the two neutral Higgs fields, and sgn(/i) is the sign of the higgsino mass parameter \i. 
With these input parameters, renormalization group equations for the CMSSM parameters 
are solved using the code Isasugra [11] to obtain the mass spectra at the weak scale. In 
the present analysis, we take A = and fi > with m t = 178 GeV (top quark pole mass) 
and mi, = 4.25 GeV (bottom quark MS mass at m&), and investigate the M 5 dependence of 
the relic density Q x h 2 in the three-dimensional CMSSM parameter space (mo, mj/2, tan/?). 

In Fig. ^ and Fig. El we show the allowed region in the (7711/2, tuq) plane consistent with 
the WMAP 2a allowed range 0.094 < £l x h 2 < 0.129. Figure [T] contains the contour plots of 
Q x h 2 for tan/3 = 50, A = 0, and /i > 0. The upper left window corresponds to the usual 
result in the standard cosmology (M 5 = 00). The dotted, dashed, solid and short-dash-long- 
dashed lines correspond to Q x h 2 = 1.0, 0.3, 0.1 and 0.05, respectively. The region among the 
bold line and the two coordinate axes is excluded by various experimental constraints (the 
lightest Higgs mass bound, b — > 57 constraint, the lightest chargino mass bound, etc.) P.ll2^ 
or the condition for successful electroweak symmetry breaking (EWSB). In particular, the 
region m <C m 1( / 2 is excluded since the lighter stau is the LSP, while the region m ^> mu2 
is excluded since successful electroweak symmetry breaking does not occur in this region 
(no- EWSB region). 

The shaded regions (in red) are allowed by the WMAP constraint. The allowed regions 
include (i) a resonance region which appears in the bulk region, (ii) a stau coannihilation 
region which appears as a narrow strip along the boundary of the stau LSP region, and (iii) 
a Higgsino-like region which appears near the no- EWSB region. In the resonance region, the 
annihilation cross section is enhanced via heavy Higgs resonances. In the stau coannihilation 
region, the lighter stau is nearly degenerated with the lightest neutralino so that the cross 
section is enhanced by coannihilation effects with the lighter stau. In the Higgsino-like 
region, the lightest neutralino is Higgsino-like, and the cross section is enhanced through 
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coannihilation effects with the second-lightest neutralino and the lighter chargino, since /x is 
relatively small in this region so that these particles are nearly degenerated with the lightest 
neutralino. 

The upper right window and the lower window in Fig. ^ are the corresponding results 
for M 5 = 4000 TeV and 2000 TeV, respectively. These figures clearly indicate that, as M 5 
decreases, the allowed regions shrink significantly. The allowed region eventually disappears 
as M 5 decreases further. 

Fig. EI is the similar result for tan/? = 30. In this case, a resonance region does not exist, 
and only narrow areas in the Higgsino-like region and the stau coannihilation region are 
allowed even in the standard case (M 5 = oo) 1 . As M 5 decreases, the allowed regions shrink 
significantly and disappear eventually as in the case of tan/5 = 50. 

Finally, we present sensitivity of the relic density to M5 in Fig. El fixing mo and mi/2 as 
well as tan/?, A and sgn(/x). The left window, where we take tan/3 = 50, m = 280 GeV 
and mi/2 = 360 GeV, corresponds to the point giving the smallest value of Q x h 2 in the 
small mass region (mo, mi/2 < 1 TeV) for tan/5 = 50. The range of Q x h 2 between the two 
dash-dotted lines satisfies the WMAP constraint. For large M5 > 10 4 TeV, the (too small) 
relic density (Q x h 2 « 0.02) in the standard case is reproduced independently of M 5 . This 
is because x t < Xd( s ) is obtained for M 5 (> 10 4 TeV). As M 5 decreases, however, the relic 
density Q x h 2 starts to be enhanced significantly, and it amounts to ~ 10 for M 5 « 100 TeV. 
It is found that, through the enhancement, an allowed region comes out for M5 in the range 
of 1000 TeV < M 5 < 1500 TeV. 

On the other hand, the right window, where we take tan/? = 30, mo = 2 TeV and mx/2 = 
490 GeV, represents the point giving the smallest value of Q x h 2 in the Higgsino-like region 
for tan (3 = 30. The M5 dependence is quite similar to the case in the left window. However, 
since the relic density Vt x h 2 rs 0.01 in the standard case is smaller than that in the left 
window, the WMAP allowed region is shifted to smaller M5 region as 600 TeV < M5 < 800 
TeV. 

In Fig. 3, we can see that fl x is proportional to about M§ 2 in the left window and M 5 ~ L6 
in the right window for a small M 5 region where neutralino abundance is enhanced. These 

1 Unfortunately, the allowed coannihilation region is too narrow to be clearly seen even in the upper right 
window. See Ref. for the figures of this region. In the other two windows, there is no allowed 
coannihilation region. 
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M5 dependences of fi x are consistent with the observation in the previous section, where 
Q x oc M 5 3 ^ 2 for S-wave (n = 0) and oc M% 3 for P-wave (n = 1). This is because in the 
left window the neutralino is Higgsino-like and its pair annihilation process is dominated 
by S-wave (n = 0). On the other hand, in the right window, the neutralino is bino-like, 
so that (for the parameters we have used in our analysis) the S-wave annihilation process 
is somewhat suppressed and the resultant M 5 dependence is expected to be in the middle 
range between M^ 1 ' 5 and M^ 3 . 

We have produced the same contour plots as Fig. [T] and El for tan/3 = 5, 10, 20 and 40 as 
well. As in Fig. ^ and El the WMAP allowed regions for each tan /3 shrink and eventually 
disappear as M 5 decreases. It is found that, almost independently of tan/3, the allowed 
region which disappears last is always the Higgsino-like region. Then, for the Higgsino-like 
region, we have plotted the same figures as Fig. El for various tan/3 = 5, 10, 20, 40 and 50. 
It turns out that the WMAP allowed range of M 5 is not so sensitive to tan/? and comes 
out around 600 TeV < M5 < 1000 TeV. Note that, since for M5 smaller than this range, 
we cannot find the WMAP allowed region, M5 > 600 TeV is the lower bound on M5 in the 
brane world cosmology based on the neutralino dark matter hypothesis, namely, the lower 
bound in order for the allowed region of the neutralino dark matter to exist. 



IV. CONCLUSIONS 



We have studied the neutralino relic density in the CMSSM in the context of the brane 
world cosmology. If the five-dimensional Planck mass is low enough, M 5 < 10 4 TeV, the 
p 2 term in the modified Friedmann equation can be effective at the decoupling time and 
the relic density can be enhanced. Therefore, we can expect that the allowed region in the 
parameter space is modified from the one in the standard cosmology. We have presented our 
numerical results and shown that the allowed region shrinks and eventually disappears as M5 
decreases. Through the numerical analysis, we have found a new lower bound on M 5 > 600 
TeV in the brane world cosmology based on the neutralino dark matter hypothesis. This 
lower bound is obtained by concerning the Higgsino-like region which disappears last. If we 
consider the resonance region or the stau coannihilation region individually for each tan /3 
fixed, we can obtain a more stringent lower bound. 
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APPENDIX A 



In this appendix, we derive an approximate formula for a decoupling temperature of 
dark matter particles in brane world cosmology and present the effect of p 2 term on the 
decoupling of the dark matter. It is shown that the decoupling temperature weakly depends 
on the transition temperature, and, thus, the resultant decoupling temperature in the brane 
world cosmology is not so much different from the one in the standard cosmology. 

Our discussion follows the strategy in We define the decoupling temperature xj, in 
terms of x, as a temperature which satisfies A(xd) = Yeq^Xcl), where A = Y — Yeq is the 
deviation of the abundance from its thermal equilibrium value. Here we examine the case 
where the annihilation process is dominated by S-wave (n = 0), for simplicity. Then, the 
condition (in the p 2 term dominated era) is expressed as 



A(x d ) ~ 

Thus, Xd is roughly written as 



— — = Y EQ ~ O.U5^-x 3 / 2 e- x 
2Ao- * g*s 



(Al) 



x d ~ In 



9 



0.145— 2Act 

9*s 







PO 


T=m 



x d ( s ) -2\nx u 



(A2) 



where 



Xd(s) = In 



0.145— 2Acx 

9*s 



(A3) 



is the decoupling temperature in the standard cosmology [9]. The second term, — 21nx t , is 
nothing but the effect of the p 2 term, and the minus sign indicates that the brane world 
effect advances the decoupling time. 



In the case that we are interested in, xm s ) < %t < m/lMeV. Using the neutralino 
decoupling temperature Xd( s ) — 30, normally used in the standard cosmology, we obtain 
7.0 < Xd < 23 for the neutralino with mass m ~ 100 GeV. The normalization used in 
Eq. (jHJ) corresponds to 

gl /A m / m \ /10 3 TeV\ 3/2 lk s 

* a V s 150 (™Giv)(^r) ' (A4) 
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FIG. 1: Contours of the neutralino relic density Q x h 2 in the (w-i^jW-o) plane for M5 = 00 
(upper left window), 4000 TeV (upper right window), and 2000 TeV (lower window) in the case of 
tan (3 = 50, A = and > 0. The dotted, dashed, solid and short-dash-long-dash lines correspond 
to Q x h 2 = 1.0, 0.3, 0.1 and 0.05, respectively. The shaded regions (in red) are allowed by the 
WMAP constraint. The region outside the bold line, including the two axes, are excluded by 
experimental constraints or the condition for successful electroweak symmetry breaking. 
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FIG. 3: n x h 2 vs. M 5 for tan/? = 50, m = 280 GeV, m 1/2 = 360 GeV (left window) and 
tan/? = 30, m = 2 TeV, m 1 / 2 = 490 GeV (right window) with A = and > 0. The range of 
£l x h? between the two dash-dotted lines satisfies the WMAP constraint. 
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